Quantum Control of a Single Qubit 
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Measurements in quantum mechanics cannot perfectly distinguish all states and necessarily dis- 
turb the measured system. We present and analyse a proposal to demonstrate fundamental limits 
on quantum control of a single qubit arising from these properties of quantum measurements. We 
consider a qubit prepared in one of two non-orthogonal states and subsequently subjected to de- 
phasing noise. The task is to use measurement and feedback control to attempt to correct the state 
of the qubit. We demonstrate that projective measurements are not optimal for this task, and that 
there exists a non-projective measurement with an optimum measurement strength which achieves 
the best trade-off between gaining information about the system and disturbing it through mea- 
surement back-action. We study the performance of a quantum control scheme that makes use of 
this weak measurement followed by feedback control, and demonstrate that it realises the optimal 
recovery from noise for this system. We contrast this approach with various classically inspired 
control schemes. 

PACS numbers: 03.67.Pp, 03.65.Ta, 03.67.-a 



I. INTRODUCTION 

Any practical quantum technology, such as quantum 
key distribution or quantum computing, must function 
robustly in the presence of noise. Many modern "classi- 
cal" technologies tolerate noise, faulty parts, etc., by re- 
lying on feedback control systems, which monitor the sys- 
tem and use this information to control its state. Given 
the ubiquity and power of feedback control for classical 
systems, it is worthwhile investigating how such control 
concepts can be applied to quantum technologies as well. 
However, strategies for quantum control must take into 
account some fundamental features of quantum mechan- 
ics, namely, restrictions on information gain, and mea- 
surement back-action. 

Classically, it is possible in principle to acquire all the 
information about the state of a system with certainty 
by using sufficiently precise measurements. That is, the 
state of a single classical system can be precisely deter- 
mined via measurement. For quantum systems, however, 
this is not always possible: if the system is prepared in 
one of several non-orthogonal states, no measurement can 
determine which preparation occurred with certainty. 

In addition, for quantum systems, monitoring comes 
at a price: any measurement that acquires information 
about a system must necessarily disturb it uncontrol- 
lably. This feature is often referred to as hack- action 

— the fundamental noise induced on a system through 
any measurement, which maintains the uncertainty rela- 
tions. This feature of quantum measurement is also dis- 
tinct from the classical situation, wherein measurements 
that do not alter the state of the system can in principle 
be performed. 

These two fundamental features of quantum systems 

— that non-orthogonal states cannot be perfectly dis- 
criminated, and that any information gain via measure- 
ment necessarily implies disturbance to the system — 



require a reevaluation of conventional methods and tech- 
niques from control theory when developing the theory 
of quantum control. 

In this paper, we investigate the use of measurement 
and feedback control of a single qubit, prepared in one of 
two non-orthogonal states and subsequently subjected to 
noise. Our main result is that, in order to optimize the 
performance of the control scheme (as quantified by the 
average fidelity of the corrected state compared to the 
initial state), one must use non-projective measurements 
with a strength that balances the trade-off between in- 
formation gain and disturbance. 

Belavkin was the first to recognise the importance of 
feedback control for quantum systems and describe a the- 
oretical framework for analysing both discrete and con- 
tinuous time models Despite this early start, it 
is only recently that the degree of control and isolation 
of quantum systems has progressed to the point that the 
experimental exploration of quantum control tasks has 
been possible H S S H S If , and the field is now un- 
dergoing rapid development (see for example 0). 

The specific control problem we are interested in here 
is the stabilization against noise of states of a single 
two level system. Similar problems have been consid- 
ered in continuous time feedback models, e.g., the sta- 
bilization of a single state of a driven and damped two- 
level atom ^1 and the maintenance of the coherence 
of a qubit undergoing decoherence p^ . Several recent 
papers have investigated state preparation and feedback 
stabilization onto eigenstates of a continuously-measured 
observable in higher-dimensional systems [la. IT^ . In 
contrast to these prior investigations, we investigate a 
feedback scheme to stabilize two non-orthogonal states 
of a two-level system. We work in a discrete-time set- 
ting, rather than continuous-time as considered in most 
prior work, which considerably simplifies the problem 
and most clearly illustrates the central concepts. Gre- 
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goratti and Werner have investigated exactly this kind of 
model of recovering the state of the system after interac- 
tion with the environment [THl fl^ in the case where it is 
possible to make measurements on the environment. In 
our setting we imagine that the environment that causes 
the initial decoherence is not available subsequently for 
the feedback protocol. Our main interest is to investi- 
gate the effects of the kind of trade-off between informa- 
tion and disturbance that is ubiquitous in quantum infor- 
mation in a concrete optimal control problem. Related 
information-disturbance trade-offs in quantum feedback 
control are discussed in ^3 • Finally, we note that imple- 
menting quantum operations on a single qubit through 
the use of measurement and feedback control as consid- 
ered here has been investigated for eavesdropping strate- 
gies in quantum cryptography^^] and for engineering 
general open-system dynamics [l9l |. 

Note that there is a fundamental difference between 
the kind of quantum control problem we are considering 
here and the related task of quantum error correction. 
(For an introduction to the latter, see |2^.) The essence 
of quantum error correction is to encode abstract quan- 
tum information into a physical quantum system and to 
choose degrees of freedom that are unaffected by the rel- 
evant noise, or upon which errors can be deterministi- 
cally corrected. However, it can be the case that one 
wishes to protect particular physical degrees of freedom 
of quantum systems and one is not free to choose an ar- 
bitrary encoding. (One such example is reference frame 
distribution via the exchange of quantum systems). The 
quantum states required for these schemes cannot be en- 
coded into quantum error correcting codes or noiseless 
subsystems _21J; protecting such systems from noise may 
therefore be an application of this kind of quantum con- 
trol. 

The paper is structured as follows. In section ^1 
we define the control task in detail; in section IIIII we 
present and determine the performance of control strate- 
gies based on "classical" concepts. Section IIVI intro- 
duces our quantum strategy, investigating the use of weak 
quantum measurements, and analyses its performance 
against the strategies of section lllll We also demonstrate 
that our quantum control scheme is optimal for the task 
at hand. In section we discuss the implications of our 
result and their relevance to other problems. 



II. A SIMPLE CONTROL TASK 

The aim of this paper is to explore the key issues we 
will confront when applying concepts from control theory 
to finite-dimensional quantum systems. In order to facil- 
itate the analysis and to be able to concentrate on the 
key departures from classical control, we will chose a very 
simple quantum system and noise model. The emphasis 
is not towards a practical task, but as an illustrative ex- 
ample. 

Consider the following operational task: a qubit pre- 
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FIG. 1: Schematic of a quantum control procedure. A qubit, 
subjected to dephasing noise, is subsequently measured and 
corrected based on the results of this measurement. The out- 
put state p is compared with the input state to charac- 
terise how well the scheme performs. 

pared in one of two non-orthogonal states \'ipi) or \'ip2) 
(with overlap ('0i|'02) = cos 6* for < 6* < 7r/2) is trans- 
mitted along a noisy quantum channel. Without knowing 
which state was transmitted, we will attempt to "correct" 
the system, i.e., undo the effect of the noise, through the 
use of a control scheme based on measurement and feed- 
back; see Fig. ^ 

The noise model that we will consider is dephasing 
noise. Let {|0), |1)} be a basis for the qubit Hilbert space, 
and the Pauli operator Z is the unitary operator defined 
by Z|0) — |0), Z\l) — — Dephasing noise is char- 
acterized as follows: with probability p a phase-flip Z 
is applied to the system, and with probability 1 — p the 
system is unaltered. The noise is thus described by a 
quantum operation j23|, i.e., a completely-positive trace- 
preserving (CPTP) map £p, that acts on a single-qubit 
density matrix p as 

£p{p)^piZpZ) + {l-p)p. (1) 

We will consider the noisy channel to be fully charac- 
terized, meaning that p is known and without loss of 
generality in the range < p < 0.5. 

We will choose the two initial states to be oriented in 
such a way that their distinguishability, as measured by 
their trace distance, is maintained under the action of the 
noise. It is straightforward to show that this condition is 
satisfied by the states 

|i/;i) = cosf|-f)-|-sinf|-), (2) 
|^2) = cosfH-)-sinf|-), (3) 

where |±) = (|0) ± |l))/\/2. 

Consider the Bloch sphere defined by states |0) and |1) 
as the poles on the z-axis. The two states and \ip2) 
lie in the x—z plane and straddle the equator of the Bloch 
sphere by angles ±6; see Fig. 12 On this Bloch sphere, the 
dephasing noise acting on these states has the effect of 
decreasing the x-component of their Bloch vectors. The 
trace distance between these two states, given by the Eu- 
clidean distance between their Bloch vectors, is invariant 
under this dephasing noise. 
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FIG. 2: Bloch sphere representation of the initial states, and 
the states after the noise. The noise shortens the Bloch vec- 
tors along the a;-axis. We have used the notation £{il') as a 
shorthand for f (I?/)) 

We now consider whether there exists a control proce- 
dure C (some "black box" ) that can correct the state of 
this system and counteract the noise, at least to some de- 
gree, independent of which input state was prepared. To 
quantify the performance of any such procedure, we will 
use the average fidelity to compare the noiseless input 
states with the corrected output states pi. Assum- 
ing an equal probability for sending either state {ipi) or 
\ip2)^ the figure of merit is 

Fc = ^Fi\^,),p,) + ^F{\^2),P2) 

= i(V'l|pi|Vl) + ^(V'2|p2|V'2), (4) 

where the fidelity between a pure state \ip) and a mixed 
state p is defined as F{\tp),p) = {ij;\p\ip). The fidelity 
F ranges from to 1 and is a measure of how much 
two states overlap each other (a fidelity of means the 
states are orthogonal, whereas a fidelity of 1 means the 
states are identical). It has the following simple op- 
erational meaning when the input state is pure: the 
fidelity F{\ipi),pi) is the probability that the state pi 
will yield outcome {ipi) from the projective measurement 

Thus, the aim is to find a control operation, described 
by a CPTP map C independent of the choice of initial 
state, such that the corrected states 

p,=C[fp(|^.)(^.|)], (5) 

for i = 1, 2 are close to the original states as quantified by 
the average fidelity. We consider control operations that 
consist of two steps: a measurement on the quantum sys- 
tem, followed by a feedback operation that is conditioned 
on the measurement result, as shown in Fig. ^ 

III. CLASSICAL CONTROL 

In this section, we introduce two types of control 
schemes for this task, both of which are based on classi- 



cal concepts, and we calculate the performance of these 
schemes based on the average fidelity. In Sec. II VI we will 
introduce a quantum control scheme that outperforms 
both of these classical schemes. 



Classical Strategy A: Discriminate and Reprepare 

For the control of classical systems, it is always advan- 
tageous to acquire as much information about the sys- 
tem as possible in order to implement the best feedback 
scheme. In line with this principle, a possible control 
strategy would be to perform a measurement on the sys- 
tem which attempts to discriminate between the input 
states, and then to reprepare the system in some state 
based on the measurement result. 

We first characterize all possible discriminate-and- 
reprepare schemes; such schemes are associated with en- 
tanglement breaking trace preserving (EBTP) maps p3 . 
|2^ . as follows. Any discrimination step is described by 
a generalized measurement, (or positive operator-valued 
measure (POVM)) ^] yielding a classical probabihty 
distribution. The generalized measurement is described 
by the operators {Pa} with Pa > and Pa = I- 
The resulting map on the quantum system is called a 
quantum- classical map QC |24| . given by 

QC{p) = J2T^T[pPa]\ea){ea\, (6) 

a 

where {|ea)} is an orthonormal basis. The reprepare 
step, in which the quantum system is re-prepared based 
on the classical measurement outcome, is described by a 
classical-quantum map CQ (24l |. given by 

CQ{p)=J2T^T[p\et){et\]Qb, (7) 

b 

where {Qb} are density matrices. 

The concatenation {CQoQC){p) leads to a map of the 
form 

B{p)^Y.^r[pP,]Q,. (8) 

b 

This map is an entanglement breaking channel. The 
name arises because the output system is unentangled 
with any other system, regardless of its input state. In 
fact it is straightforward to see from [2^ |2^ that all 
EBTP maps can be realised by some discriminate-and- 
reprepare scheme. Thus these EBTP maps formalize our 
notion of discriminate-and-reprepare strategies. 

The measurement for discriminating two (possibly 
mixed) preparations given by Helstrom 25] is optimal 
in terms of maximizing the average probability of a suc- 
cess. For our choice of states, Helstrom's measurement is 
a projective measurement onto the basis {|0), |1)}, which 
successfully discriminates the states and \^2) with 
probability PHcistiom = ^(1 + sin0). Note that because 
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of the particular choice of dephasing noise, this success 
probabihty is independent of the noise strength p. 

We now present and analyse two possible discriminate- 
and-reprepare strategies, both of which are based on Hel- 
strom's measurement. 



Discriminate and Reprepare Scheme 1: 

With the outcome of Helstrom's measurement, one 
strategy is to reprepare the qubit in either state \ipi) or 
I ■02) based on this measurement outcome. This scheme 
yields an average fidelity of 

Fum = 1-5 (sin^ 6 - sin^ 0) . (9) 

Such a replacement ignores the fact that the discrimina- 
tion step can fail, with probability 1 — PHeistrom, in which 
case a prepared state lipi) would be reprepared as \4>2) 
(or vice versa). 



Discriminate and Reprepare Scheme 2: 

We can consider other strategies that reprepare differ- 
ent states so as to reduce the effect of the aforementioned 
error. In particular, we now demonstrate that the follow- 
ing pair of states maximises the average fidelity: 

l*±) = + VsT^ll) , (10) 

where 7 = \/ sin'* 6 + cos^ 9. Note that this replacement 
is also independent of p. Here, is prepared if the 

measurement outcome corresponds to \ipi), and l^*-) 
is prepared otherwise. In this strategy, the reprepared 
states are slightly biased towards the alternate state to 
that suggested by the measurement (smaller 9) — in a 
sense hedging our bet. As a proof of the superiority of 
this scheme over the former, the fidelity 



i^DR2 - - + -Vcos2 0-Hsin''0, (11) 

satisfies Fdr2 > ^dri for all 9. Both Fdri and i^DR2 are 
presented in Fig.|3fa). 

This second discriminate-and-reprepare scheme is in 
fact the optimal discriminate-and-reprepare scheme, in 
that it achieves the highest average fidelity 

2 

maxFg = maxi^(VdS[£:p(|V^)(V'^l)]l^z), (12) 

4=1 

where the maximization is over all EBTP maps B acting 
on a single qubit. This optimization was performed (in 
a different setting) by Fuchs and Sasaki ^26]. In the Ap- 
pendix, we provide an alternate proof of optimality using 
techniques from convex optimization. 



Classical Strategy B: Do Nothing 

Another control strategy would be to do nothing to 
correct the states. Although trivial, this strategy is of 
interest for comparison with other schemes. (There ex- 
ist schemes that perform worse than this strategy, be- 
cause of the feature of quantum systems that every mea- 
surement that acquires information will uncontrollably 
disturb the system.) This scheme does not lie within 
the set of discriminate-and-reprepare schemes described 
above (it is not described by an EBTP map) but we will 
nonetheless refer to it as "classical." 

The average fidelity of this scheme is given by 

Fn = 1 -pcos^e. (13) 

This performance is plotted in Fig. |SIb). Clearly, this 
scheme performs best for small amounts of noise (p ~ 0) 
and for input states with Bloch vectors that are near the 
z-axis (which is invariant under the dephasing noise). In 
some non-trivial regions of the (p, 9) parameter space, 
in particular in the range of low noise, this "do noth- 
ing" scheme outperforms the optimal "discriminate and 
reprepare" scheme. 

IV. QUANTUM CONTROL 

In the previous section, we presented control schemes 
based on classical concepts. However, using techniques 
that may lead to optimal control schemes for a classical 
system may not necessarily lead to optimal schemes for a 
quantum system. As we will now demonstrate, the above 
classical control strategies can be outperformed by using 
a strategy based on quantum concepts. 

We note that the two classical schemes presented in the 
previous section lie at the extreme ends of a spectrum: 
the "discriminate-and-reprepare" strategy achieved max- 
imum information gain and induced a maximum distur- 
bance, whereas the "do-nothing" strategy achieved zero 
disturbance but produced zero information gain. As 
demonstrated by Fuchs and Peres [23|, there exist an 
entire range of generalized measurements that trade off 
information gain and disturbance. A possible avenue for 
improvement in our control schemes is to tailor the mea- 
surement in such a way as to find a compromise, if one 
exists, between acquiring information about the noise but 
not disturbing the system too much result of the 
measurement. 

In the following we re-express the noise process £p in 
a way that suggests a strategy for constructing such an 
improved feedback protocol. 

A. Reexpressing the noise 

To develop an intuitive picture, we will make use of 
a preferred ensemble for the quantum operation £p de- 
scribing the noise. That is, we use a decomposition of 
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a) b) c) 




FIG. 3: The performance of the schemes, quantified by the average fidehty, as a function of the amount of noise p and the angle 
between the input states 9. a) Discriminate and reprepare scheme quantified by the average fidelity -Fdr2 of Eq. Illlll . The 
fidelity Fori of Eq. @ is shown as a solid line at p = 0.5. Both average fidelities -Fdr2 and -Fori are independent of p. b) "Do 
nothing" scheme, quantified by the average fidelity _Fn of Eq. 1131 . For this scheme, the average fidelity drops to Jn = 1/2 for 
p = 0.5 and 6 = 0. c) Quantum control scheme, quantified by the average fidelity -Fqcopt of Eq. 12711 . The range of fidelities 
plotted has been made identical in all the figures to aid comparison. 



the operation into different Kraus (error) operators than 
that given in Eq. The resulting quantum operation 
£p describing the noise, however, is equivalent. 

Consider the following quantum operation on a qubit, 
viewed on the Bloch sphere: with probability 1/2, the 
Bloch vector of the qubit is rotated by an angle +a 
about the z-axis, and with probability 1/2 it is rotated 
by —a about the z-axis. Rotations about the z-axis are 
described by the operator 

= e"'"^/^ = cos(a/2)/ - i sin(a/2)Z , (14) 

and the quantum operation is then 

= sm^{a/2){ZpZ)+cos^{a/2)p. (15) 

Thus, this quantum operation is equivalent to the de- 
phasing noise £p, with p = sin^(a/2). 

Viewing the noise operation £p with this preferred en- 
semble, it is possible to describe the noise as rotating the 
Bloch vector of the state by ±a with equal probability. 
A possible control strategy, then, would be to attempt to 
acquire information about the direction of rotation (±a) 
via an appropriate measurement, and then to correct the 
system based on this estimate. Loosely, we desire a mea- 
surement that determines whether the noise rotated the 
state one way (-l-a) or another (—a). Then, based on the 
measurement result, we apply feedback: a unitary oper- 
ation (rotation) that takes the state of the system back 
to the desired axis. 

A projective measurement, wherein the state of the 
system collapses to an eigenstate of the measurement, 
does not meet these requirements because such a mea- 
surement destroys the distinguishability of the two pos- 
sible states. Instead, we consider the use of a weak mea- 
surement, with a measurement strength chosen to bal- 
ance the competing goals of acquiring information and 
leaving the system undisturbed. We now show that such 
a strategy is possible, and that there is a non-trivial op- 
timal measurement strength for this task. 



B. Weak non-destructive measurements 

For our quantum control scheme, we will make use 
of a type of measurement that satisfies two key require- 
ments: (1) the strength of the measurement should be 
controllable, i.e., we should be able to vary the trade-off 
between information gain and disturbance (back-action); 
and (2) the measurement should be non-destructive, 
which leaving the measured system in an appropriate 
quantum state given by the desired collapse map. Such 
weak non-destructive measurements have recently been 
developed and demonstrated in single-photon quantum 
optical systems [23, ll^ . 

Using the preferred ensemble describing the noise, 
Eq. H15|l . we expect intuitively that this weak measure- 
ment should be along the y-axis of the Bloch sphere in 
order to provide information about which direction (±a) 
the system was rotated, without acquiring information 
about which initial state the system was prepared in. 
One suitable family of POVMs consists of two operators 
given by Em = Mj^Mm, for m ~ 0, 1, where M,„ are the 
measurement operators |23| 

Mo = cos{x/2)\+i)(+i\+sm{x/2)\-i){-i\ , (16) 
Ml = sin(x/2)|+*)(+*| + cos(x/2)hz)(-i| . (17) 

The strength of the measurement depends on the choice 
of the parameter x- The eigenstates of Y are |±«) = 
(|0) ± i\l))/^/2. The probabilities of obtaining the mea- 
surement results TO = 0, 1 for a qubit in the state pin are 
given by 

Prn = Tr[ErnPin] , (18) 

and the resulting state of the qubit immediately after the 
measurement is 
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Consider the following two limits. If x = 7r/2 the 
two measurement operators are the same and are pro- 
portional to the identity. As a result the outcome prob- 
abilities are independent of the state and the state of 
the signal is unaltered by the measurement. If x = 0, 
a projective measurement on the signal is induced: the 
signal state is projected onto the state when 
the measurement result is (1). For < x < 7r/2, the re- 
sulting measurement on the signal is non-projective but 
non-trivial. 

It is illustrative to view the effect of this measurement 
on the noisy input states on the Bloch sphere. In Fig.^Ja) 
we can see that the effect of the noise is to shorten the 
length of the Bloch vector of the qubit state (making it 
less pure) while increasing the angle between the Bloch 
vector and the x-y plane from 6 to 6', where 6' > 9 . 
When the measurement is made, three things happen, as 
can be seen in Fig.^fb): 1) the Bloch vector is lengthened 
(the state becomes more pure); 2) the angle 9' decreases 
to some lesser angle 9"; and 3) the state is rotated about 
the z-axis one way or the other depending on the result 
of the measurement. The first two effects work towards 
our advantage (purifying the state while decreasing 9'); 
the third effect we attempt to correct using feedback. 

We will now describe how to implement this measure- 
ment using a projective measurement on an ancillary me- 
ter qubit and an entangling gate between the original sig- 
nal qubit and the meter. The strength of the measure- 
ment can be controlled by varying the level of entangle- 
ment between the two qubits, which can be implemented 
by initiating the meter in the state |0) and subsequently 
applying a rotation (as shown in figure [SJa)), where 



|0) 



V - _ ( cos(x/2) - sin(x/2)\ 

" " Uin(x/2) cos(x/2) i 



(20) 



The parameter x ranges from to 7r/2 and characterises 
the strength of the measurement, with equivalent to a 
projective measurement and 7r/2 equivalent to no mea- 
surement. 

The entangling gate consists of a Xjl rotation on the 
signal state, followed by a CNOT gate with the signal state 
as the control and the meter state as the target, followed 
by a X^iL on the signal state, where 



-i(t>X/2 



cos(0/2) 
—i sin((/)/2) 



-'isin((?!)/2)\ 
cos(0/2) ) 



(21) 



and where the Pauli matrix X is given by X\Q) — |1) 
and X\l) — |0). The rotations X±]l are used to ensure 
that the resulting weak measurement on the signal qubit 
is performed in the basis. The entangling 

gate then correlates (to a degree which depends on x) 
the {|-f«), \—i)} basis of the signal qubit to the {|0), |1)} 
basis of the meter qubit. 

Finally the meter qubit is measured in the basis 
{|0),|1)}, yielding a result or 1. This measurement 
on the meter induces a measurement on the signal that is 
precisely equal to the generalized measurement described 
by the measurement operators Mm of Eq. H16|) . 
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FIG. 4: Bloch sphere representation of the effect of a weak 
measurement on the system. The transformations shown here 
correspond to having obtained the measurement result "0" 
(for the result "1", the behaviour would be a reflection in the 
x-z plane.) a) The two initial states \^iii,2) are mapped to 
Pi, 2 by the noise; b) a weak measurement is performed with 
< X < f i c) a strong projective measurement (x = 0) is 
performed projecting either state into \—i)- While no mea- 
surement will not yield any information about the system, a 
strong measurement will maximally disturb the system. A 
weak measurement will gain some information while also lim- 
iting the disturbance on the system. 



C. Feedback control 

Once a weak measurement has been performed, a cor- 
rection based on the measurement result is performed on 
the quantum system: the feedback control. We choose 
the correction to be a unitary rotation about the z-axis, 
Z±ri where 



-irjZ/2 



c- 





-vn/2 



(22) 



with the aim to bring the Bloch vector of the qubit back 
onto the xz-plane. The angle of rotation is chosen to 
be ±77, depending on the measurement result (+ry cor- 
responding to the measurement result 0, and —77 to the 
measurement result 1). It is possible to choose 77 so that 
the system state is returned to the a;z-plane for all val- 
ues of p, 9 and x and for both measurement outcomes by 
choosing 



tan?7 



1 



(1 — 2p) cos tan X ' 



(23) 



with T] in the range < < 7r/2. This angle rj can be cal- 
culated because the dephasing noise has been previously 
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Quantum control scheme 
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FIG. 5: a) Circuit diagram of the weak measurement scheme. 
The input signal state pin is entangled to the meter state us- 
ing the CNOT gate. The X±:^ rotations ensure that the weak 
measurement of the signal state is made in the desired basis 
|— i)}. The strength of the measurement is set using 
the rotation Y^. The meter state is measured in the compu- 
tational basis, resulting in a classical signal (0 or 1) to be fed 
forward to the correction stage of the control scheme, 
b) Circuit diagram of the control scheme. A weak measure- 
ment is made on the input state and, based on the measure- 
ment results, the signal state will be rotated by Z,, (^-i?) 
conditional on the result of the weak measurement being 
(1). 



characterised (i.e., p is known). 

The resulting weak measurement followed by feedback 
is thus described by a quantum operation (a CPTP map) 
Cqc acting on a single qubit, given by 

Cqc(p) = {Z+r,Mo)p{Z+r,MQ)^ + (Z_„Afi)p(Z_,,Mi)t , 

(24) 

where the measurement operators are given by 

Eqs. (cani). 

In summary, the quantum control scheme operates by 
performing a weak measurement of the system and then 
correcting it based on the results of the measurement, 
as in Fig. |5j3) . The weak measurement is made by en- 
tangling an ancillary meter state with the signal state 
using an entangling unitary operation, then performing 
a projective measurement of the meter state. The level 
of entanglement depends on the input state of the meter, 
which is controlled by a rotation; this level of entan- 
glement in turn determines the strength of the measure- 
ment. After measurement of the meter, the signal state is 
altered due to the measurement back-action. To correct 
for this back-action, a rotation about the z-axis is applied 
to the state, returning it back to the xz-plane. To char- 
acterise how well the scheme works, we now investigate 




Strength 



1 . 



FIG. 6: (a) Fidelity of the quantum correction procedure 
with measurement strength (1 — 2x/ti) for a representative 
noise value (p = 0.145) and angle (6' = 0.715). The measure- 
ment strength ranges from a value of (corresponding to no 
measurement), through to a value of 1 (corresponding to a 
projective measurement). There exists an optimum measure- 
ment strength at which we balance the amount of information 
gained with the amount of back-action noise introduced. Also 
plotted for comparison are (b) the optimal "discriminate-and- 
reprepare" scheme and (c) the "do nothing" scheme for the 
same parameter values. 



the average fidelity. 



D. Performance 



The performance of this quantum control scheme, 
quantified by the average fidelity Q , is 



Fqc — 



1 + sin-^ sin X + cos 6* W 1 - (1 - r^) sin x 



(25) 

where Tx = (1 — 2p) cos is the x component of the Bloch 
vector describing the system after the noise. 

We can see that Fqc is a function of the amount of 
noise p, the angle between the initial states ^, and the 
measurement strength x- The dependence of this fidelity 
on the measurement strength, for fixed p and 0, is illus- 
trated in Fig. El For each value of p and 6, there is an 
optimum measurement strength Xopt which maximizes 
the average fidelity (|25|l . This optimum measurement 
strength is found to be non-trivial except for the limiting 
cases of p = or 6* = 0, 7r/2, and is given by 



Xopt(p, 0) = sin" 



y (l-r2)2cos2 6l + (l-r2)sin^6'' 

(26) 

as a function of the amount of noise p and the angle 
between the initial states 0. 

Substituting Xopt for x in Eq. f^ . we get the follow- 
ing expression for the optimum fidelity: 



QCopt 




COS^ 



(27) 



8 




0.4 - 0.005 



0.5l ^ ^ : . ^ . ^ 1 I— lO 

n/2 n/i 

e 



FIG. 7: A contour plot of the difference, as a function of the 
amount of noise p and the angle between the initial states 6, 
between the average fidelities of the quantum control scheme 
and the best classical scheme. The quantum control scheme 
performs significantly better for moderate values of p and 8 
(0.05 < p < 0.3 and 0.3 < 6* < 1). The maximum value 
Fdif = 0.026 occurs at p = 0.115 and 9 = 0.715. 

Fig. I^c) plots the quantum control fidelity as a func- 
tion of the input state (characterised by the angle 9) and 
the amount of noise (characterised by p) . 

We note that -Fqcopt = 1 for three limiting cases. If 
p = 0, there is no noise and so the state is not perturbed, 
resulting in unit fidelity for all values of 9 given by sim- 
ply "doing nothing" (zero measurement strength and no 
feedback). When 6 = tt/2, the states are orthogonal and 
point along the z axis. The noise does not affect these 
states, again resulting in unit fidelity for all values of p 
with a "do nothing" scheme. When 9 = the two states 
are equal and point along the a:-axis. The control scheme 
reprepares this state after the noise by making a projec- 
tive measurement x = to obtain either | +i) or | —i) and 
rotating back to the xz-plane (rj = tt/2). This results in 
a fidelity of 1 for all values of p. 



E. Comparison with Classical Schemes 

We now compare the quantum control scheme with 
classical schemes presented in Sec. IIIII Specifically, we 
compare the quantum scheme with the best of the classi- 
cal schemes at every point in the parameter space (p, 9) , 
i.e., we observe the difference in the average fidelities 

^dif = -Fqcopt - max(^bR2, ^n) , (28) 

where JbR2 and Fn are given by Eqs. lfTT|l and ((T^ . 
respectively. Fig. [7| reveals that Fdn is always positive, 
and thus the quantum control scheme always outperforms 
the best of the classical strategies. 



F. Optimality 

We now prove that our quantum control scheme is 
optimal, in that it yields the maximum average fidelity 
of all possible quantum operations (CPTP maps). Our 
proof makes use of techniques from convex optimization 
(specifically, those of 'sS'l) but is presented without re- 
quiring any background in this subject. In the Appendix, 
we provide a more detailed construction of the proof. 

Consider the following optimization problem: deter- 
mine the maximum average fidelity 

2 

Fopt =maxFc = max i ^ (tA^IC [£'p(|i/',)(-0,|)] iV'i) , 
1=1 

(29) 

where the maximization is now over all CPTP maps C 
acting on a single qubit. 

Recall that any CPTP map C acting on operators on 
a Hilbert space 7i is in one-to-one correspondence with 
a density operator Tc onH($Ti. with 

C{g) = Tri„ [{g^ ® I)Tc] , (30) 

and is subject to the constraint Trout [T"c] = I'm, where 
'in' denotes the first subsystem and 'out' denotes the 
second |^ [s^l- With this isomorphism, the av- 
erage fidelity Fc for the control scheme C is given by 
Fc = Tr[i?Tc], where 

2 

i?= IV'.OV'^I- (31) 

i=l 

Thus, the optimization problem (|29|1 can be rewritten as 
maximize Tr [RTq] 

subject to Tc > (32) 
Trout [Tc] = /in ■ 

We now wish to prove that the maximum value of 
Tr [i?Tc] subject to these constraints is given by -Fqcopt 
of Eq. lf77|l . 

We note that, for any single-qubit operator M satisfy- 
ing M ^ I — R> 0, we obtain the inequality 

Tr [M] - Tr [i?Tc] = Tr [(M ® /)Tc] - Tr [RTc] 
= Tr [(M R)Tc] 
>0, (33) 

where the first line follows from the constraint 
Trout [Tc] = /in, and the inequality follows from the fact 
that (M (g) / - i?) > and Tc > 0, and thus the trace 
of their product is non-negative. This inequality demon- 
strates that the value Tr [M] for any matrix M that sat- 
isfies the constraint (M (g) / — i?) > provides an upper 
bound on the solution of our optimization problem l|32|) . 
Consider the matrix M = bo{I + rxX), where 

^o = | + J^cos2^+|^, (34) 
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and r^. = (1 — 2p) cos9 as before. It is straightforward to 
verify that the matrix h^I ® I + r^h^X ^ I — R> 0, and 
hence the value Tr [M] — 2bo provides an upper bound 
on the average fidehty of any control scheme. Because 
2bo precisely equals the fidelity of our proposed quantum 
control scheme, given by Eq. (|27(l . this scheme necessar- 
ily gives an optimal solution to the original problem (|32|l . 
We refer the reader to the Appendix for a more construc- 
tive proof of this result. 



V. DISCUSSION AND CONCLUSIONS 

We have shown how two key characteristics of quan- 
tum physics - that non-orthogonal states cannot be per- 
fectly discriminated, and that any information gain via 
measurement necessarily implies disturbance to the sys- 
tem - imply that classical strategies for control must 
be modified or abandoned when dealing with quantum 
systems. By making use of more general measurements 
available in quantum mechanics, we can design quantum 
control strategies that outperform schemes based on clas- 
sical concepts. In particular, we have presented a task 
for which the optimal scheme relies on a non-trivial mea- 
surement strength, one that balances a tradeoff between 
information gain and disturbance. 

In constructing our quantum control scheme for the 
particular task presented here, we made use of sev- 
eral intuitive guides. First, we used a preferred (and 
non-standard) ensemble of the dephasing noise operator 
(Eq. (|15|l 'l. which allowed us to view the noise as "kick- 
ing" the state of the qubit in one direction or the other 
on the Bloch sphere. We then made use of a weak mea- 
surement in a basis that, loosely, attempted to acquire 
information about the direction of this kick without ac- 
quiring information about the choice of preparation of the 
system. It is remarkable (and perhaps simply lucky) that 
these intuitive guides lead to a quantum control scheme 
that was optimal for the task. It is interesting to consider 
whether such intuition can be applied to quantum control 
schemes in general, and if this intuition can be formalized 
into rules for developing optimal control schemes. 

While our scheme is indeed optimal for the task pre- 
sented, it is not guaranteed to be unique; in fact, there 
are other decompositions of the same CPTP map into 
different measurements and feedback procedures 33] . In 
general, it is possible that an entire class of CPTP maps 
may yield the optimal performance. Also, the intuitive 
guides discussed above for our quantum control scheme 
— such as that the measurement essentially gains infor- 
mation only about the noise and not the choice of initial 
state — may not apply to other optimal schemes. 

In connection to this, we note that a similar feedback 
control scheme was investigated by Niu and Griffiths 
for optimal eavesdropping in a B92 quantum cryptogra- 
phy protocol see also In their scheme, the aim 
of the weak measurement was to maximize the informa- 
tion gain about which of two non-orthogonal states was 



transmitted for a given amount of disturbance; in con- 
trast, our weak measurement was designed to acquire no 
information about the choice of non-orthogonal states. 
Despite these opposing aims, the obvious similarity be- 
tween these our scheme and that of Niu and Grifhth 
warrants further investigation, particularly since we note 
that optimal feedback protocols exist based on different 
choices of measurement. 

Finally, we note that the key element to our quan- 
tum control scheme — weak QND measurements on a 
qubit, and feedback onto a qubit based on measure- 
ment results — have both been demonstrated in recent 
single-photo n q uantum optics experiments. Specifically, 
Pryde et al [23 have demonstrated weak QND measure- 
ments of a single photonic qubit, and have explicitly 
varied the measurement strength over the full parame- 
ter range. Also, Pittman et al |35j | have demonstrated 
feedback on the polarization of a single photon based 
on the measurement of the polarization of another pho- 
ton entangled with the first; this feedback was used for 
the purposes of quantum error correction, and is essen- 
tially identical to the feedback required for our quantum 
control scheme. Because these core essential elements 
have already been demonstrated experimentally, we ex- 
pect that a demonstration of our quantum control scheme 
is possible in the near future. 
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APPENDIX: OPTIMIZATION PROOFS 

In Sec. lIIII and llVI the proposed classical and quantum 
control schemes were shown to be optimal among the set 
of EBTP and CPTP maps, respectively. Here, we provide 
constructive proofs of these results in further detail. 

1. Weak duality 

Consider the following optimization: 
maximize Tr [-fo^] 

subject to Z > (A.l) 

where the matrices Fq, Fi and the vector c are spe- 
cific to the problem, and Z is the variable over which 
the optimization is performed. We say that any Z sat- 
isfying the constraints of the problem is feasible. An 
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optimization problem of this form is known as a semi- 
definite program (SDP), a class of convex optimization 
problems ;36] . Each problem of the form (|A.1(I has a La- 
grange dual optimization problem that arises from using 
the method of Lagrange multipliers and has the form ,36.] 



minimize c x 

subject to — Fo + XiFi > 



(A.2) 



where now the vector x is the variable to be optimized. 

In many cases, such as the optimization problems in- 
vestigated here, the dual problem is straightforward to 
solve, or else efficient numerical solutions are known that 
solve the primal and dual problems together. Further- 
more, the dual problem allows us to bound the optimum 
of the original problem and this fact can be used to prove 
the optimality of solutions as follows. 

Let t) = c^x be the value of the objective function to be 
minimized in ljA.2|) for an arbitrary feasible x. Similarly, 
let p = Tr [FqZ] for an arbitrary feasible Z and let p* be 
the optimum of our original problem ljA.l|) . 

We now demonstrate that, if one can find a feasible 
point to (|A.1I) yielding p and a feasible point to l|A.2|) 
yielding such that = p, then p = p*, that is, the point 
Z yielding p is optimal. 

Consider the difference 



t) — p = c X 



-Tr[FoZ]-Tr[(^ 



F,X,,-Fn]Z 



(A.3) 



where we have used the linearity of the trace and c^x = 
J2i '^i^i — Si Tr [FiZxi] . As the trace is over the product 
of two positive semi-definite matrices, it has to be non- 
negative. That is to say that 



P < P* < 5. 



(A.4) 



Clearly, if there is a p such that p = J) for some 0, then 
p = p*. 

2. Dual optimization for quantum control 



As demonstrated in Sec. IIVFI obtaining the maxi- 
mum average fidelity can be expressed as the optimiza- 
tion problem (|32|l . For this problem (as for the classical 
problem which we address in the next section) the dual 
optimization proves to be straightforward to solve ana- 
lytically and the results above can then be used to show 
optimality of the control scheme given by Eq. H24|) . 

We make use of some symmetry arguments to sim- 
plify the problem. This optimization problem has cer- 
tain symmetry properties under the action of the group 
of transformations generated by the rotation T (X (g) 
X)T{X (g) Xy and the transpose T T-^. Specifically, 
the objective function is invariant under the action of 
this group since Ty[R{X X)T{X X)] = Tr[i?T] and 
Tt[RT^] = Tr[i?T], because {X ® X)R{X ® X) = R 
and i?"^ = i?, respectively. In addition, the constraints 



are covariant under the action of the group. Since 
conjugation with a unitary and transposition preserve 
eigenvalues, {X ® X)T{X ® X) > and > if 
T > 0. To see that the equality constraints are covariant 
note that Trout [f c] — hn is equivalent to the condition 
Tr [{M ® /out)T] = TrAf for aU hermitian M. If T obeys 
the partial trace constraint we have 



Tr[(M®/out)(^®^)T(X(; 

= Ti[{XMX < 



'X)] 

5/out)T] = TrM, (A.5) 



and 



Tr[(M (g Jout)T^] = Tr[(M^ ® /out)T] = TrM , (A.6) 

so both {X (g) X)T{X (g) X) and do also. So both the 
objective function and the feasible set of (|32|l are invari- 
ant under the action of the group. As a result there will 
be an invariant point T*^^ = {X ®X)T^^^{X ®X) = T*^ 
that achieves the optimum p* [s^. We do not need to op- 
timize over the full set of T but may restrict our attention 
to the set of invariant Tinv Gatermann and Parrilo j37j 
have investigated such invariant SDP's in detail. 

The dual of our optimization problem l|32|) has the 
form [ai 



minimize TrM 
subject to M ® I ■ 



R>0 



(A.7) 



Notice that (as is generally the case) this semidefinite 
program is invariant under the same group of transforma- 
tions as the original problem, under which M XMX 
and M — > . For the dual problem we may likewise 
restrict attention to Mjnv = bgl + bxX that are invariant 
under the action of the group. This gives a simpler dual 
optimization 



minimize 26o 
subject to bol ( 



> I + b^X (g) / - i? > . 



(A., 



where bo and b^ are the new variables. This problem is 
simple enough to solve analytically; the solution is 



1 



/cos2 6l 



sm 



l-r? 



(A.9) 



and bx = r^bQ (with r^ = (1 — 2p)cos6). This may be 
checked by verifying that the matrix b^I ® I ~\-bxX ® I — R 
is indeed positive semi-definite, hence 26o is a valid dual 
feasible value. Because 2&o reproduces the fidelity of our 
proposed scheme, given by Eq. this guess neces- 

sarily gives an optimal solution to the original problem 
(123. 



3. Dual optimization for classical control 

The same approach is used to solve the problem H12|l . 
We start by mapping the set of trace-preserving entan- 
glement breaking qubit channels to bipartite states Tg. 
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For these channels Tg is positive, has partial trace equal 
to the identity, and is also separable (2^]. Because Tg 
is an (unnormalised) state of two qubits, the separabil- 
ity condition is equivalent to the positivity of the partial 
transpose 38]. We will denote the partial transpose of 
the operator Tg on the subsystem Tiout by T^°"' . Thus 
we may rephrase the optimization problem (|12|l in the 
form 



maximize Tr [i?Tg] 
subject to T/3 > , T; 

Trout Tr = In 



> 



(A.IO) 



Note that the condition of positivity of the partial trans- 
pose guarantees that Tg corresponds to an EBTP map. 

The new problem has the same symmetries as the full 
optimization H32(l with one addition. Notice that = 
R so the objective function of both problems is invariant 
under partial transpose. In our new problem the point 
Tg°"' is feasible if Tg is feasible, so the feasible set is also 
invariant under the partial transpose. (Note that since 
partial transpose does not preserve positivity this is not 
true of the problem H32I) '). Because of this symmetry 
we may restrict our attention to Tjnv for which Tj^"^" = 
Tinv Since the partial transpose sends A^Y^ —A (g) 
Y where A is any Hermitian matrix, we can conclude 
that Tt[{A (g) F)Tinv] =0. It is sufficient to check this 
condition for the full set of Pauli matrices /, X, Y, Z so 
the requirement of invariance under the partial transpose 
constitutes four new constraints. Notice however that 
> is now redundant since we are 
jj^^ = T"^'"^. So we can replace the 
problem ljA.10|) with 



the condition T^°^ 
requiring that T ' 



maximize Tr [i?Tg] 
subject to Tg > 

Trout Tg = -fin 

Tr(A(gy)Tg = VAe {I,X,Y,Z}. 

(A.ll) 

Positivity of the partial transpose and hence the separa- 
bility of Tg is now guaranteed by the positivity of Tg 
and the additional equality constraints. 
The dual of the problem IjA.lip is 



minimize TrM 
subject to / (g) M -h iV ( 



Y-R>0 



(A.12) 



This semidefinite program still has symmetries corre- 
sponding to the rotation X ^ X and the transpose (but 
not under the partial transpose.) These two symme- 
tries lead to the transformations N —XNX and 
N —N'^ respectively. The only invariant choices of 
N are proportional to Y. As before we may restrict at- 
tention to Minv — aol + axX that are invariant under 
the action of the group and iVinv = ayY . This gives a 
simpler dual optimization 



minimize 



2ao 

subject to flo/ I 



/ -f- a^X g) I - 



.Y®Y 



where ao, ax and a„ are the new variables. 



i? > 0, 
(A.13) 
This problem 



should be compared to the analogous dual optimization 
in the quantum case ljA.8|l . Again, this problem can be 
solved analytically, yielding the solution 



ao 



- + - \/ cos2 B + sin* B , 
4 4 

rx cos^ 
4 4 v/cos2 e + sin'* e 
cos 9 sin^ 9 



4 V 



(A.14) 
(A.15) 

(A.16) 



cos^ 



sin 



Again, one can check that agl g) / -I- UxX g) / -I- UyY g) 
y — i? is positive semidefinite with these choices, which 
ensures that the objective function 2ao is indeed a dual 
feasible value. The proof of optimality follows as before in 
the quantum case by: (i) observing that 2ao reproduces 
the fidelity i^DR,2 of Eq. Hll|l and (ii) applying the weak 
duality argument. 

We note that the optimization techniques presented 
here may be useful when applied to more general prob- 
lems presented in Fuchs and Sasaki |2^. However, when 
the map in question does not act on qubits, there are sig- 
nificant complications in characterizing the EBTP maps 
because the PPT condition is no longer sufficient. 
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